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ABSTRACT 

The model considered is a d = 2 layered random Ising system on a square lattice 
with nearest neighbours interaction. It is assumed that all the vertical couplings 
are equal and take the positive value J while the horizontal couplings are quenched 
random variables which are equal in the same row but can take the two possible values 
J and J — K in different rows. The exact solution is obtained in the limit case K — > oo 
for any distribution of the horizontal couplings. The model which corresponds to this 
limit can be seen as an ordinary Ising system where the spins of some rows, chosen at 
random, are frozen in an antiferromagnetic order. No phase transition is found if the 
horizontal couplings are independent random variables while for correlated disorder 
one finds a low temperature phase with some glassy properties. 
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Ising spin glasses have been solved exactly in their mean field version [1,2] while as 
far as I know no exact solutions are available at finite dimensionality d > 2. Indeed, 
in presence of disorder even a d = 1 system with magnetic field is a very complicated 
problem [3] and compact exact solutions can be found only in special cases [4]. This 
is frustrating since it is not always clear if the qualitative results of a mean field 
approximation are shared by the finite dimension model. In this paper I am far from 
answering this problem, nevertheless I exactly solve a class of d = 2 layered random 
Ising systems which in some conditions have a low temperature phase. The nature 
of this non-ferromagnetic low temperature phase is still unclear to me, but there are 
some indications that it shares some of the properties of a glassy phase. 

The models I consider are defined as follows: the interaction is effective only 
between nearest neighbours on a square lattice; all the vertical couplings are equal 
while the horizontal couplings are quenched variables which are equal in the same row 
but can take different values in different rows; the vertical couplings take the positive 
value J while the horizontal couplings J, can take the two possible values J and J — K 
with K — > oo. These models have frustration since the product of the signs of the 
coupling around a plaquette can be negative. Different distributions of the horizontal 
couplings correspond to different models of the class; in the simplest case the Jj are 
independent random variables and take the value J with probability 1—p and J — K 
with probability p. 

Layered Ising models of this type have been first considered by B.M. McCoy 
and T.T. Wu [5,6] for the non frustrated case which is used for studying the effect 
of quenched randomness on the ferro-para transition. These authors deal with the 
determinant which occurs in the Pfaffian approach and while they do not provide an 
explicit exact solution of the problem they are able to show that the free energy has an 
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infinitely differentiable singularity at the transition. Layered models with frustration 
have been studied by R. Shankar and Ganpathy Murthy [7], not only their topic but 
also their approach is the same of this work since they deal with the row to row 
transfer matrices. They do not find out an exact solution, nevertheless they map the 
problem into a collection of d = 1 random field Ising systems from which they can 
extract a lot of informations. In particular they provide evidence for the existence of 
a low temperature phase. 

Let me now state more precisely the problem. Assume that N = LM is the number 
of spins, L is the number of rows and M the number of columns, the hamiltonian can 
be written as 

h n = (Jaija i+1 j + Ji<Ti,j<Ti,j+i) (1) 

ij 

where the Jj are the horizontal couplings whose value only depend on the row i and 
not on the column j. One can write Ji = J — rjiK where the quenched variables rji 
can take the value and 1 according to a given distribution. In the independent case 
rji = with probability 1 — p and rji = 1 with probability p. The partition function is 

Z N = ^2exp(^2/3(Jaija i+lij + Jaija iJ+1 - rjiK{l + aija iJ+1 ))} (2) 

{o-} ij 

where the constant term £V rjiK has been added to the hamiltonian in order to avoid 
divergences in the K — > oo limit. After having defined Y = J(3 and performed the 
limit K — > oo one obtains 



Zn = yi n ex p{ rc7 ^+i,j + ^1,3^,3+1} - 1 + at, ^ i,j+i 



(3) 



{<t} ij 

The terms in parenthesis equal 1 when rji = and (1 — o~ijo~ij + i)/2 when rji = 1. 
Notice that in this second case the antiferromagnetic order between neighbour spins on 
the row is imposed, in fact, if o^j and o~i,j+i have the same sign they give a vanishing 
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contribution to the partition function. It is now clear that (3) defines a class of Ising 
model with both vertical and horizontal couplings equal to J and with the spins of 
some rows frozen in an antiferromagnetic order. The frustration comes out from the 
fact that the tendency to the ferromagnetic alinement due to the positive couplings is 
in competition with the tendency to the antiferromagnetic alinement induced by the 
frozen spins on the unfrozen ones. A similar problem, where the spin are randomly 
frozen in a random direction has been solved in d = 1 in [4], and studied in d = 2 at 
zero temperature in [8]. 

The advantage of considering layered disorder is that one can apply a standard 
diagonalization method [7,9], and reduce the problem to the evaluation of the trace 
of products of random matrices. Following the same steps of [7,9] one easily finds the 
free energy 

/ = ~ log(2 shih2r) 7(<Z, T)dq (4) 

where 

1 L 

1 (q,T) = lim -logTrTTT^r) (5) 

L—>oo L 

i=l 

The 2x2 matrices Ti(q 7 T) can be written as the product Ti(q 7 T) = Ei -T(q, T) where 
the 

*-(J A) 

are random and equal the identity when rji = and the up projector r + = (1 + r^)/2 
when r\i = 1. The matrix T(q, T), on the contrary, is constant and reads 

T(q,T) = exp{-rV 3 } exp{2r*(r 3 cos q + n sing)} exp{-rY 3 } (7) 

where T\ and r 3 are Pauli matrices and r* = — | log(tanhT). 
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The trace of a product of random matrices is easily accessible via computer 
simulation but it cannot be, in general, exactly computed. In the present case, 
nevertheless, following a similar method as in [4] , it is possible to find out the compact 
analytical result. Consider a given realization of the quenched variables r\i and look 
at the product of matrices in (5). Since T^q, T) = ■ T(q, T) that product reduces 
to a product of matrices T(q,T) and up projectors r + . The first and the second 
t + will be separated by h matrices T{q, T), the second and the third by h matrices 
T(q, r), and so on. The l n are random variable which can take the values 1,2,.... 
whose distribution can be easily found out once the distribution of the rji is given. 
The order number n goes from 1 to rif = L/l, in fact, one must have Xm=i In = L so 
that Xm=i ln/nf = 1 = L/rif. With the help of these considerations one can rewrite 
(5) as 



where [T(g,r) z ™]n is the up left entry of T(q,T) ln . If P(l) is the probability that 
two successive rows of infinitely negative couplings are separated by / rows of finite 
positive couplings than 1 = Yl'ili IP (I) and 



7 (g,r)= lim ilogTT[T(g,r) z »] n = lim \ Vlog[T( ? ,r) z "] 
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(8) 






In order to find the explicit form for (9) it is convenient to write 



T(q, T) = exp{e(r3 cos + T\ sin 4>) } 



(10) 



where 



cosh e = 



cosh 2 2r 



(11) 



sinh 2r 



— cosq 



and 



COS = 



cosh 2r(cos q — sinh 2r) 



(12) 



(sin 2 q + cosh 2 2r(cosg - sinh2r) 2 )^ 
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Using (10) it is immediate to obtain 

[T(g,r)']n = cosh(/e) + cos0sinh(Ze) (13) 

Finally: 

/ = — -log(2sinh2r) ~J^P(l) / log(cosh(/e)+cos0sinh(/e))dg (14) 

2T 2nTl Jo 

where e and <p are given in (11) and (12). 

The probability P(l) for the simplest choice of independent rji is P{1) = p(l — p) 1-1 
and I = 1/p. In this case one can prove that the system has no phase transition, except 
for p — where it trivially reduces to the ordinary Ising model. In Fig. 1 it is shown 
the specific heat C in correspondence of different values of p; one can notice that the 
logarithmic divergence is smoothed showing the absence of transition. Nevertheless, 
the model is frustrated and its zero temperature properties are not completely trivial. 
One can compute the T — energy /o end entropy sq and finds 

f = -2J(l-p) 2 (15) 

s = Jp\l-p)log(^±^j (16) 

so is not vanishing for p ^ 0,1 showing an exponential degeneration of the ground 
state due to the frustration of the model. 

Since the transition disappears for p 7^ the role of p reminds that of a 
magnetic field which also suppresses the transition. The analogue of the spontaneous 
magnetization is obtained in the limit p — > as 




This quantity is continuous while its derivative ^ is not as shown in Fig. 2 where 
one can see a logarithmic divergence at T c (the Onsager critical temperature). 
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The circumstance that a phase transition can be found only at p = suggests 
to look more carefully at the model around this value. If one chooses p = a/L one 
has a vanishing p in the thermodynamic limit and the free energy is the same of 
that of the standard Ising model. Nevertheless, one has a random finite number of 
frozen rows. This number is Poisson distributed with intensity a and it is different 
for different realizations of the disorder (no self-averaging). The distance between 
two given frozen rows is also a random number of order of L and it also varies from 
a realization to another. The final result is that the frozen rows separate a random 
number of regions of random size of order N whose magnetization at T < T c is ±m(T) 
independently one from the other (m(T) is the Onsager spontaneous magnetization 
at temperature T). As a consequence of this fact, the whole system can be in all the 
states corresponding to all the possible combinations of magnetization of each region. 
In conclusion, one has the same free energy of a standard Ising system but a number 
of pure states each of them corresponding to a different local magnetization. The 
situation is completely analogous to that studied in [10] for a diluted d = 1 model at 
zero temperature. Following the same line of [10] it is easy to compute the overlap 
probability, in particular for large a one has 

P(,),-i= exp{-|!} (18) 

where t = m(T) 2 /a. (18) implies that the overlap between two different pure states 
vanishes in the limit a — > oo , nevertheless it should be noticed that the self-overlap 
Qmax = m(T) 2 , being independent on a, remains finite. 

To summarize: for p ^ there is no phase transition, while for p = a/L one has a 
glassy like phase which comes out from an artificial construction which maintain the 
same free energy of the Ising model. It is straightforward, at this point, to look at an 
intermediate situation, where the number of frozen rows is of order of L but they can 



be much more separated than in the independent case. This task can be accomplished 
with the choice P(l) = a /I s (a is the normalization constant) which replaces the 
exponential distribution P(l) = p(l — p) l ~ x of the independent case. By substituting 
this expression in (14) one can easily compute the free energy and look at the eventual 
divergences. In spite of the fact that the free energy is now different from that of the 
standard Ising model one still finds a phase transition at the Onsager temperature. 
Nevertheless this phase transition does not correspond to a divergence in the specific 
heat, but in its derivative which is plotted in Fig. 3. I have not been able to 
quantitatively characterize the low phase temperature with an order parameter. In 
this phase, in fact, while the spontaneous magnetization vanishes, both the overlap 
and a parameter connected with the antiferromagnetic order seems to differ from zero. 
I hope that some light on this point will come out from future research. 
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Figure Captions 



Fig. 1 Specific heat C as function of the temperature T. The dotted line corresponds to 
the Ising model (p = 0), the full line to p = 0.1 and the dashed line to p = 0.2 



Fig. 2 Temperature derivative of /' = 



dp 



as function of T 



p=0 



Fig. 3 Temperature derivative 4£ of the specific heat for the P(l) = a/l 3 model. 
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